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Abstract 

The simultaneous QCD analysis of the xF^ structure functions measured in 
deep-inelastic scattering by several collaborations is done up to 3-loop order 
of QCD. The x dependence of the higher-twist contribution is evaluated 
and turns out to be in a qualitative agreement with the results of "old" 
CCFR data analysis and with renormalon approach predictions. The Gross- 
Llewellyn Smith sum rule and its higher-twist corrections are evaluated. 
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1. The experimental data of the CCFR collaboration (we'll call them "old") ob- 
tained at Fermilab Tevatron for the xF^ structure functions of deep-inelastic scat- 
tering of neutrinos and antineutrinos on an iron target provide an important means of 
accurate comparison of QCD with experiment. However, in view of revision of "old" 
data announced dy CCFR collaboration P] the question arises: what can we say about 
the comparison of the QCD predictions on dependence of the xF^^XjQ"^) structure 
function (SF) based on the data of neutrino DIS experiments different from those of 
CCFR? 

In the present note, a combined fit of the experimental data of the CDHS [Q, SCAT 
g, BEBC-WA59 §, BEBC-Gargameile [§ and JINR-IHEP collaborations for the 
xF^ structure functions is done in order to determine the x dependence of the SF, higher 
twist (HT) contribution and the value of the scale parameter A;^. 



2. We'll use, for the QCD analysis, the Jacobi polynomial expansion method pro- 
posed in 1^. It was developed in [§]-[|l4| and applied for the 3-loop order of pQCD to 
fit F2 and XF3 data [|T§ 0. 

The - evolution of the moments M^'^^^ [N, Q^) is given by the well known per- 
turbative QCD |T^, |T^ formula: 



pQCD 



H^{Ql,Q')Mf''''{N,Ql), N = 2,3,... (1) 



The factor Hm (Qq, Q"^) contains all next- and the next-to-next-to-leading order 



QCD corrections □ and is constructed in accordance with |TJ] based on theoretical results 
of @. 

The expression (|I|) provides an input for reconstruction of the SF by the Jacobi 
polynomial method. Following the method [|10|, 0, we can write the structure function 
xF^ in the form: 

^ max TT' 

xFf'^'^x, Q') = x-il - xf E Q-/{x) E cf\m?'"' (j + 2, Q') , (2) 

n=0 i=0 

where 9"^(x) is a set of Jacobi polynomials and c^{a,l3) are coefficients of the series 
of G"'^(x) in powers of x: 



(3) 



j=0 

The unknown coefficients M^lN, Qq) in (|I]) could be parametrised as Mellin moments 
of some function: 



Mi 



pQCD 



(N, Ql) 



dxx^-^Ax\l - xy{l + 7x), N = 2, 3, ... 

Jo 



(4) 



^For reviews and references on higher order QCD results see[[l8[ 
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To extract the HT contribution, the nonsinglet SF is parameterized as follows: 



xF^ix, Q'') = Q'') + h{x)/Q'' 



(5) 



where the dependence of the first term in the r.h.s is determined by perturbative 
QCD. Constants h{x.i) (one per x-bin) parameterize the HT x dependence. We put 
= 0.03, 0.05, 0.08, 0.15, 0.25, 0.35, 0.45, 0.50, 0.55, 0.65, 0.80 for i = 1, 2. ..11. 
The HT contribution for F2 was determined in pO|. The values of constants h{xi) as 
well as the parameters A, b, c, 7 and scale parameter A are determined by fitting the 
combined set of data of 192 experimental points of xF^ in a wide kinematic region: 
0.5 GeV^ < < 196 GeV^ and 0.03 < x < 0.80 and Ql = 10 GeV^. We have 
put the number of flavors to equal 4. In accordance with the result of concerning the 
disagreement of their data with perturbative QCD at small x, a cut x > 0.35 was used 
for CDHS data. The TMC are taken into account to the order o^M^^^J Q'^) . 

The nuclear effect of the relativistic Fermi motion is estimated /^from below by the 



ratio R 



D/N 



F^/F^ obtained in the covariant approach in light-cone variables |21 





LO 


NLO 


NNLO 


Xd.f. 


312/176 


316/176 


312/176 


A 


6.68 ± 0.38 


6.92 ± 1.43 


7.11 ± 0.38 


b 


0.760 ± 0.027 


0.768 ± 0.072 


0.778 ± 0.027 





4.03 ± 0.07 


3.97 ± 0.17 


3.82 ± 0.07 


7 


0.675 ± 0.156 


0.452 ± 0.624 


0.189 ± 0.128 


Am5 [MeV] 


191 ± 46 


159 ± 39 


163 ± 31 


Xi 


h{xi) [GeV^' 




0.03 


0.086 ± 0.087 


0.090 ± 0.091 


0.067 ± 0.085 


0.05 


0.001 ± 0.028 


0.022 ± 0.032 


0.093 ± 0.047 


0.08 


-0.127 ± 0.123 


-0.094 ± 0.126 


-0.011 ± 0.131 


0.15 


-0.286 ± 0.046 


-0.230 ± 0.050 


-0.200 ± 0.050 


0.25 


-0.401 ± 0.058 


-0.334 ± 0.056 


-0.327 ± 0.054 


0.35 


-0.284 ± 0.073 


-0.220 ± 0.068 


-0.178 ± 0.062 


0.45 


-0.436 ± 0.093 


-0.366 ± 0.090 


-0.403 ± 0.083 


0.50 


0.005 ± 0.079 


0.047 ± 0.077 


0.036 ± 0.074 


0.55 


-0.243 ± 0.069 


-0.200 ± 0.068 


-0.242 ± 0.064 


0.65 


0.176 ± 0.063 


0.202 d= 0.072 


0.154 ± 0.060 


0.80 


0.020 ± 0.037 


0.024 ± 0.039 


-0.012 ± 0.039 



Table I. Results of the 1-, 2- {NMax = 10) and 3- order {NMax = 8) 
QCD fit (with TMC) of the combined xF^ SF data for / = 4, > 
0.5GeV^ with the corresponding statistical errors, normalization 
coefficients and values of the HT contribution h{xi). 

3. Results of the fit are presented in Table 1 and Figures 1-3. The theoretical 
prediction for h{x) from |2^ is presented Figure 3. 

The experimental values of xF^ for each collaboration were multiplied by the nor- 
malization factors C"^°" which were considered as free parameters. Their values are not 
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sensitive to the order of pQCD in use and was found to be equal to: C^^^'" ^^^^ = 
0.92 ± 0.03, C^^^^ = 1.06 ± 0.03, C^^^^-^^^^ = 1.02 ± 0.05 and C'^^^^-^"^^'- = 
0.97 ± 0.04. The value of C^^^s ^ ^ g^ed. 

The obtained value of Ajjo is larger than that given by a similar analysis of CCFR 



data ||T5[ ^jjg = 134 ± 57 MeV but exhibits relatively small statistical errors. Results 
of the NLO and NNLO fit give the constant of strong interaction a^'^'^(M|) = 0.105 ± 
0.004 and ag^^'^{M'^) = 0.107 ± 0.003 in agreement, within the errors, with usual 
DIS results |]2^. Additional uncertainties to the value of as(M|) due to extrapolation 
of the dependence of the SF with four flavors (f=4) in a wide kinematic interval 
0.5 GeV^ < < 196 GeV^ were found to be 0.001 in |6l . 

13] at = 10 GeV^ 
-- 2.60 ± 0.23 in 



The value of the perturbative part of the GLS sum rule 



1 



dx 



estimated by using results of Table 1 is equal to /q 
agreement with results of the "old" CCFR data analysis |^ 

The shape of h(x) is in qualitative agreement with theoretical predictions of the 
dispersion method of the renormalon approach ( for reviews and references see p3| ) 
and with results of the QCD analysis of "old" CCFR data presented in |]T5[. They 
obviously differ from the precise values of h(x) for singlet ST F2 presented in |2^ . 

Based on the results of Table 1, one can estimate the value of the first moment of 
h(x) which contributes to the GLS sum rule 



27 : hi 



The obtained values: 



h{o = -0.42±0.27 H , h^^^ = -0.29±0.28 and h^^^^ = -0.26±0.27 are in agreement 
with theoretical predictions of ^ h = -0.29 ±0.14 and §^ h = -0.47 ±0.04 as well 
as with the recent result of pl| . 



4. In conclusion it should be stressed that combined fit provides still a more precise 
determination of A-j^ and h{xi) in comparison to the analysis of "old" CCFR data 



T5| , while the shape of the SF ruled by parameters A, b, c and 7 is determined less 
accurate. The most discrepancy with the "old" CCFR data analysis takes place for the 
HT contribution to the GLS sum rule and for the HT x dependence at large x. 

For a more precise determination of the HT contribution to SF, the role of the nu- 
clear effect should be clarified and a more realistic approximation for Rp^^^ = F^^ / 
is needed. We also did not take into account the threshold effects on Q"^ evolution of SF 
due to heavy quarks P2| which is necessary owing to a wide kinematic region of data un- 



der consideration and could be realized based on the mass- dependent MOM-scheme [26 
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•^Hereafter present the value of h(x) in [GeF^]. 
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Figure captions. 

Fig.l. Higher-twist contributions from LO fit. 
Fig. 2. Higher-twist contributions from NLO fit. 

Fig.3. Higher-twist contributions from NNLO fit and the theoretical prediction for 



h{x) from ||22 . 
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